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A critical look at V^s determinations from hadronic r decay data 

Kim Maltman, Math and Statistics, York University, 4700 Keele St., Toronto, ON CANADA M3J 1P3* 

A critical review of hadronic r decay data based determinations of \V US \ is given, focussing on the impact of 
the slow convergence of the integrated D — 2 OPE series for the conventional flavor-breaking sum rule deter- 
mination and the potential role of as-yet-unmeasured multiparticle contributions to the strange spectral distri- 
bution. Additional information obtainable from analyses of inclusive strange decay data alone, and from mixed 
electroproduction-r sum rules with much reduced OPE uncertainties, is also discussed. Self-consistency tests are 
shown to favor determinations which reduce somewhat discrepancies with 3-family unitarity expectations. 



1. Introduction 

Kfz and T[K — > nv\/T\it — > pv] based deter- 
minations of \V US \ [T] are in excellent agreement 
with expectations based on \ V u d\ [2] and 3-family 
unitarity. In contrast, determinations from fla- 
vor breaking (FB) sum rules involving hadronic 
t decay data [3141516] yield results ~ 3<r low, even 
after accounting for recent HFAG strange exclu- 
sive branching fraction (BF) updates [6]. 

For a kinematic-singularity-free correlator, n, 
with spectral function, p(s), the finite energy sum 
rule (FESR) relation on which the FB r determi- 
nation is based takes the form 

/■so i r 

/ w(s)p(s) ds = <z> w(s)H(s) ds 

Jo 2tti /| s | =ao 

(1) 

a result valid for any sq and any analytic w(s). 
\V US \ is obtained by applying Eq. (JTJ) to the FB 



difference An, 



n 



(J) 

VjA-M 



n 



(0+1) 
V+A-ud 



n 



(0+1) 
V+A;us 



where 



(s) are the spin J = 0, 1 components of the 
flavor ij, vector (V) or axial vector (A) current 
two-point functions, and (0 + 1) denotes the sum 
of J — and 1 components. For sufficiently large 
,s , the RHS of Eq. (TTJ) is evaluated using the OPE 
representation, [An r ] OPB , while, for s < m\, 
the LHS is obtainable from inclusive hadronic 
t decay distributions. Explicitly, the spectral 
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functions, 



Pv/A 



are related to the differential 



distributions, dRv/A;ij/ds, of the normalized fla- 
vor ij V or A current induced decay widths, 

Rv/A-i] = rb"~ -» v T hadronsvyyt;ij (7)]/ r [ T ~ -> 
v T e~v e (j)], by [7] 



dR 



V/A-ij 



„EW 



ds 



-WL{yr)Pv/A;iji S ) 



(2) 



with y T = a/ml, w r{v) = (1 - J/) 2 (l + 2y), 
wl{u) = 2y(l — y) 2 , Vij the flavor ij CKM matrix 
element, and, with Sew & short-distance elec- 
troweak correction [8], d? w = \2-k 2 S ew / m2 ■ 

The choice of FESRs involving the J = + 1 
combination, An T , rather than that correspond- 
ing to the linear combination of spectral func- 
tions appearing in Eq. @, is predicated on the 
extremely bad behavior of the integrated J = 0, 
D = 2 OPE series for scales kinematically accessi- 
ble in r decay [9] . Fortunately, the dominant J = 
spectral contributions are the accurately known, 
non-chirally-suppressed tt and K pole terms. The 
remaining, continuum contributions, are (i) dou- 
bly chirally suppressed (Pv+A-ij K l( m i -f m j) 2 ]), 
(ii) numerically negligible for ij = ud, and (iii) 
determinable phenomenologically for ij = us J = 
via dispersive [10] and sum rule pj] analyses 
of the strange scalar and pseudoscalar channels. 
The J = contributions can thus be subtracted 
bin by bin from dRv+A-,ij/ds, allowing one to de- 
termine Pv+A-ij(. a ) ano - cons truct the re- weighted 
J = + 1 spectral integrals, Ry +A . i j{so), and FB 
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differences, SR^ +A (s ), denned by 



R V+A:ij( S 0) 

6RV +A (s ) -- 



R V+A;ud( S o) _ Ry+A;us( S o) 

\V ud \ 2 \V US \ 2 
dsw(s)Ap T (s) . (3) 



EW 



Taking | V u d | and any parameters in the OPE rep- 
resentation of 8Ry +A from other sources, Eq. ([1]) 
then yields [3] 



IK. I = 



\ 



n 



V+A;us 



(so) 



|K,d| a 



X r>w,OPE, \ 

oR v+A (s ) 



(4) 



For weights used previously in the literature [3] 
1415] . the OPE contribution to the denominator 
is at the few-to-several-% level of the ud spec- 
tral integral term, making modest accuracy for 
fiRy+A E ( s o) sufficient for a high accuracy de- 
termination of \V US \- The smallness of the FB 
OPE corrections is illustrated by the results of the 
,s = m 2 ., w T FESR where, with updated HFAG r 
BFs [BJ, supplemented by Standard Model (SM) 
K^2 expectations for the r — > Kv T BF, Bk, the 
value \V US \ = 0.2108(19) obtained ignoring FB 
OPE corrections differs from that obtained using 
one of the possible evaluations of 5Ry T _^ A PE (m 2 ), 
\V US \ = 0.2174(22) [6J, by only - 3%. 

The task of quantifying the uncertainty in the 
OPE estimate for 5Ry T +A (s ) and, from this, the 
theoretical error on \V US \, is complicated by the 
rather slow convergence, at the correlator level, of 
the series for the leading dimension D = 2 OPE 
contribution [An r ]^^. To four loops [12] 



ope _ 3 m s (Q 2 ) 
D = 2 ~ 2^2 Q2 

19.93a 2 + 208.75a 3 + d 4 a 4 



[AH T (Q 2 )] 



(5) 



with a = a s (Q 2 )/n, and a s (Q 2 ) and m s (Q 2 ) the 
running coupling and strange quark mass in the 
MS scheme q Since a(m 2 ) ~ 0.1, convergence 
at the spacelike point on the contour \s\ = sq 



2 In what follows, we employ the estimate d± ~ 2378 1121 
for the as-yet-undetermined 5-loop coefficient d^. 



is marginal at best. With such slow convergence, 
conventional error estimates may significantly un- 
derestimate the D = 2 truncation uncertainty. 

Fortunately, the FESR framework allows for 
internal self-consistency checks. Assuming both 
the data and OPE error estimates are reliable, 
the \Vus\ obtained from Eq. (j4|) should be inde- 
pendent of So and w(s). On the OPE side, three 
common methods can be employed for evaluating 
the integrated D = 2 contribution: the contour 
improved (CIPT) prescription, used with either 
(i) the truncated expression for [AIT r ] I)=2 itself, 
or, (ii) after partial integration, with the corre- 
spondingly truncated Adler function, and (iii) the 
truncated fixed-order (FOPT) prescription (for 
which the correlator and Adler function results 
are identical). While the CIPT prescription cor- 
responds to the conventional procedure of choos- 
ing a local scale, arguments in favor of the FOPT 
alternative also exist |13j . At a given truncation 
order, the three prescriptions differ only by con- 
tributions of yet higher order. If the estimated 
truncation errors are reliable, the results for \V US \ 
obtained using the different prescriptions should 
thus also agree within these errors. 

2. \V US \ from FB hadronic r decay FESRs 

The results in this section are obtained using 
the updated December 2010 HFAG t BFs [6], 
supplemented by SM expectations for Bk- 
SM 7r M 2 expectations are similarly employed for 
the non-strange 7r pole contributions. Assuming 
lepton universality, the constrained HFAG deter- 
mination of the electron BF, B e = 0.17852(27), 
implies Ry+A-ud = 3.467(9) and R v +A;us = 
0.1623(28). The change in Rv+A-.ud necessitates 
a small rescaling of the publicly available 2005 
ALEPH ud distributions [l4j , from which the con- 
tinuum contributions to Pv/A-ud{s) are obtained. 
For the continuum part of pv+A-,us, the BaBar 
and Belle analyses of the inclusive strange distri- 
bution are not yet complete. The completed 1999 
ALEPH distribution [15] corresponds to exclusive 
strange BFs with significantly larger errors, and, 
in many cases, significantly different central val- 
ues, than those obtained by the B factory ex- 
periments [B], As an interim measure, we fol- 



3 



low the strategy of Ref. [T^] , "partially updating" 
Pv+A-us ( s ) to reflect the new values of the BFs 
through a mode-by-mode rescaling of the 1999 
ALEPH distribution. This proceedure, though 
not ideal, has been tested using the BaBar dis- 
tribution data for r — > K~ 7r + ir~ v T [17] (a mode 
having a particularly large BF change) , and found 
to yield spectral integral contributions in surpris- 
ingly good agreement with those of the actual 
BaBar data. The rescaling method, however, has 
not been checked for other modes, and is likely 
to be less reliable for the X37r, KAtt, ■ ■ ■ contri- 
butions, where the ALEPH distribution was not 
measured, but estimated using Monte Carlo, and 
has a purely phase space shape. 

OPE input is as in the last of Refs. [4], except 
for the update a" /=3 (m?) = 0.3181(57), which 
reflects the new world average, o£ f ~ (M%) = 
0.1184(7) [IE]. 

We begin with the conventional w T determina- 
tion. For so = "v, the ud and us spectral inte- 
grals are determined by the corresponding inclu- 
sive BFs. Improvements to the various exclusive 
strange BFs [6] then translate into an improved 
determination of the us spectral integral, even 
without completion of the remeasurement of the 
inclusive strange distribution. This is not the case 
for other so and/or other w(s). The D = 2 trun- 
cation uncertainty in this case has been estimated 
using standard last-term-retained©residual-scale- 
dcpcndcncc methods (known to work well for per- 
turbative series displaying good convergence be- 
havior). The resulting combined theoretical un- 
certainty on \V US \ is then 0.0005 [5]. To test for 
the independence of \V US \ on sq and w(s), one 
must go beyond w T and/or sq = m^- We perform 
these tests using the interim updated pv+A-us(s)- 

In Fig. [TJ results for \V US \ as a function of sq 
are shown for each of the three prescriptions for 
the w r -wcightcd D = 2 OPE series. The two 
CIPT-bascd prescriptions give similar results, but 
show poor so-stability. The FOPT prescription 
yields significantly improved, though not perfect, 
so-stability. For all s the FOPT and CIPT re- 
sults differ by significantly more than the 0.0005 
total theoretical error estimate mentioned above. 

The observed s -instability of the FB w T FESR 



- - w , CIPT+correlator 
_ w*, CIPT+Adler 

— w ,FOPT 

v 




Figure 1. \V US \ vs. so from the FB uv-wcighted 
FESR for the three different prescriptions for han- 
dling the integrated D = 2 OPE series. 



results could arise from deficiencies in any of the 
contributions appearing in Eq. @. Problems 
with the ud spectral integrals seem unlikely, but 
can be tested for using FESRs for the us V+A 
channel alone. A problem with missing higher 
multiplicity us spectral strength could also be ex- 
posed by such an analysis, assuming the OPE rep- 
resentation used to be reliable. Specifically, for So 
large enough to include some portion of the region 
of missing strength, |V^ S | would be too low, while 
for so below the region of the missing strength, 
\V US \ should rise back to its true value. The us 
V+A FESRs are considered in the next section. 

The most obvious candidate for a source of so- 
instability in the OPE contribution is the poten- 
tially slowly converging integrated D = 2 series. 
For so = m^, the w T -weighted integrated D = 2 
series behaves as l+0.29+0.10-0.04-(0.20) + - • • 
for the CIPT+Adler function case, 1 + 0.15 + 
0.02 - 0.12 - (0.29) + • • • for the CIPT+correlator 
case, and 1 + 0.40 + 0.26 + 0.15 + (0.08) + • • ■ 
for the FOPT case, where the terms in paren- 
theses are the estimated 5-loop, 0{a 4 s ) contribu- 
tions. For the two CIPT cases, though there is 
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cancellation on the contour at intermediate or- 
ders, this cancellation does not persist to higher 
orders, and the integrated series appears rather 
badly behaved. The behavior is better for the 
FOPT case, which was also found to yield im- 
proved So-stability. Worth noting is the fact that 
5R 



,OPE 



V+A 



mr T ) for the better behaved FOPT pre- 
scription is a factor of ~ 2 larger than that for ei- 
ther of the two CIPT prescriptions. The - 100% 
increase between CIPT and FOPT increases \ V MS \ 
by ~ 0.0020 establishing the insufficiently con- 
servative nature of the 0.0005 estimate for the 
combined theoretical error. With the difference 
between the FOPT result and the average of the 
two CIPT results as a new truncation uncertainty 
estimate, this source dominates the theoretical er- 
ror. Taking the FOPT result (favored by the sq- 
stability criterion) as the new central value, the 
result for the FB w T determination becomes 



|K s |=0.2193(3) ud (19) us (19) t7l 



(0) 



~ 2.3a below 3-family unitarity expectations. 

FB r-based FESRs can also be constructed us- 
ing other w(s). For the CIPT prescription, if the 
so-histability of the w T FESR results from pre- 
mature truncation of the D = 2 series, conver- 
gence can be improved by choosing weights which 
emphasize contributions from the region of the 
complex s = —Q 2 plane away from the spacclikc 
point, where |a s (Q 2 )| is smaller and the conver- 
gence, at the correlator level, of the D — 2 scries 
better. Three weights (denoted u>io(y), wio{y), 
and u>2o(y), with y = s/sq) having this, as well 
as other desirable properties, were constructed in 
Ref. [19]. The results for \V US \ from the result- 
ing FB FESRs, using the CIPT+correlator pre- 
scription, are shown in Fig. [2j The w T results 
for both the FOPT and CIPT+correlator pre- 
scriptions are shown for comparison. Improved 
so-stability compared to the w T CIPT results is 
observed^. The \V US \ values are compatible with 



3 FOPT-based results for the new weights arc not shown, 
though they also display significantly improved so- 
stability. The reason is that, having been constructed to 
improve the convergence of the integrated CIPT series, 
there is no reason to expect the new weights to also im- 
prove the convergence of the integrated FOPT series, and 
indeed they do not. For example, the FOPT version of the 



the result of Eq. ^ within the theoretical er- 
rors quoted previously for that case. The central 
result for wio (the weight showing the best sq- 
stability) is \V US \ = 0.2188. The experimental 
error, which, absent a re-measurement of the in- 
clusive strange distribution, has to be based on 
the 1999 ALEPH us covariances, is 0.0033. 



— w , FOPT 

. . . w , CIPT+corr 



w 10 , CIPT+corr 



w 10 , CIPT+corr 




s [GeV ] 



Figure 2. \V US \ vs. sq, obtained using the 
CIPT+correlator D = 2 OPE prescription for 
the FB wio, wio and W20 FESRs, and, for com- 
parison, using FOPT and CIPT+correlator treat- 
ments for the w T FESR. 



3. The us V+A and mixed t-EM FESRs 

FESR analyses of the us, J = + 1, V+A cor- 
relator combination provide an alternate deter- 
mination of \V US \, one that provides some cross- 
checks on the FB FESR results. Two new ingredi- 
ents enter the OPE side of these FESRs: D = 
contributions, determined by a s , and a D = 4 
gluon condensate contribution absent from the 



wiio-weighted D = 2 series behaves as 1 + 0.51 + 0.43 + 
0.39 + (0.41) -I , in contrast to the CIPT+correlator ver- 
sion, which behaves as 1 + 0.24 + 0.19 + 0.15 + (0.11) H . 
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FB difference AI1 T . It is known that, with the 
5-loop D = expansion as input |20j . pinched 
FESR analyses of the ud, J = + 1 V, A and 
V+A correlator combinations [2T] yield values 
of a s in excellent agreement with current high- 
precison lattice determinations [22j and the new 
world average [18] . A strong anti-correlation be- 
tween a s and the gluon condenstate |21| implies 
(a s G 2 /n) = 0.012 GeV 4 for the central a s input 
noted above. Results for \V US \ as a function of sq 
obtained from the iu T -weighted FESR using each 
of three prescriptions for the integrated D = 2 
series are presented in Fig. [3J The compatibil- 
ity of the results for each prescription with those 
of the corresponding FB w T FESR is excellent. 
However, the so-dcpcndcnce of |Ki S | for the two 
CIPT prescriptions is clearly not compatible with 
the assumption that the D = 2 OPE representa- 
tion is reliable and the instability in the FB w T 
FESR results for \ V US | due solely to missing higher 
multiplicity us spectral strength. As for the FB 
FESR case, the FOPT D = 2 treatment produces 
improved, though not perfect, so-stability. 



0.225 




s [GeV ] 



Figure 3. \V US \ vs. sq for the u> r -based us V+A 
FESR, using the three different prescriptions for 
handling the integrated D = 2 OPE series. 



With problems in the FB r FESRs, to at least 
some extent, due to slow D = 2 OPE conver- 
gence, FESRs with D = 2 OPE contributions 
reduced at the correlator level are highly desir- 
able. With this in mind, combinations of nj^ 1 ]^, 

Hv+A-us an d tne EM correlator, Hem, (whose 
spectral function, pem, is determined by the bare 
e + e~ — > hadrons cross-sections) can be con- 
structed having zero D = and vanishing O(ctg) 
D = 2 OPE contributions [23]. The unique such 
FB combination with the same normalization for 



n 



(o+i) 

V+A;us 



as An r is 



An A/ = m EM - 6n£$ + An T 

The D = 2 OPE series for A±T M (<2 2 ) is 



3 m s 
2^2 Q2 



4.38a 2 + 44.9a 3 



(7) 



which has not only (by construction) a vanish- 
ing O(a s ) coefficient, but also higher order coef- 
ficients significantly smaller than those for An r , 
c.f. Eq. ([S]). The resulting strong suppression 
of the D = 2 series is accompanied by a fortu- 
itous suppression of D = 4 contributions. Explic- 
itly, with 64 = (m s ss — niiii), the D = 4 series 
for An T and ATIm have the form ^ J2k 
with {cq,Ci,c 2 ) = (-2,-2,-26/3) for AIL, and 
(0,8/3,59/3) for AII M Integrated D > 4 

contributions, which are expected to be some- 
what enhanced [23], can be fitted to data using 
their stronger so dependences. 

A complication for the AHm FESRs is the 
fact that EM and r results for the 2ir and 
47r components of the / = 1 V spectral func- 
tion are not compatible at the level of expected 
isospin-breaking corrections [25] . We nonetheless 
present, as an illustration, the results for \V US \ as 
a function of sq for the w T , W2(y) = (1 — j/) 2 and 
W3{y) = 1 - |y + \y z An A/ FESRs, obtained as- 
suming the r data to be correct in both cases. Be- 
cause of the current experimental complications, 
no attempt has been made to fit D > 4 OPE 
contributions. Instead, D = 6 contributions have 
been estimated using the vacuum saturation ap- 
proximation and D = 8 contributions neglected. 
The main point is to illustrate the significantly 
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improved ,s -stability produced by the strong sup- 
pression of the integrated D = 2,4 OPE contri- 



butions. The Sq 



result in this case is 



|V US | = 0.2222(20) T (28)js M (9) 

where only the experimental errors have been dis- 
played. While the very small theoretical errors for 
the mixed t-EM FESRs make them promising for 
the future, further progress requires a resolution 
of the experimental 2tt and Air problems. 



0.23 
0.225 



0.22 



0.215 



- T-EM, w 

- t-EM, w 2 

- T-EM, w 3 

. . ud-us w , CIPT+Adler 




s [GeV ] 



where R TtS = Rv+A-us and R t ,v+a = Rv+A,ud 
are the i? e -normalized inclusive strange and 
non-strange branching fractions, defined above 
Eq. @. The theoretical contribution, 5R T jh, is a 
sum of J = (longitudinal) and J = + 1 (L+T) 
contributions, denoted SR T \ L and 5R T jh\ L+T , re- 
spectively, in Rcf. [5]. The former, which must 
be determined phenomcnologically because of the 
problematic behavior of the longitudinal OPE 
representation, is dominated by the accurately 
known K and n pole contributions. Adding the 
current best assessment of the small continuum 
contributions yields the result 

5R T \ L = 0.1544 ± 0.0037 (11) 

quoted in Rcf. [5]. 

SR T ,th\ L+T is evaluated using the OPE. The 
OPE representation is a sum of contributions, 
SR OPE;D> of dimensions D = 2,4, 6 • • • . D > 6 
contributions are typically neglected in the liter- 



ature. SR 



L+T 

OPE\D=2 



is proportional to m s (2 GeV) 



and, for S EW = 1.0201 [8], m T = 1.77677 ± 
0.00015 GeV 0, m s (2 GeV) = 94 McV, and 
the current nj = 5 world average a s (Af|) = 
0.1184±0.007 [IE] (equivalent to the n } = 3 result 
a s {m T ) = 0.3181 ± 0.0057) has the values given 
in Table [TJ as a function of the truncation order, 
tit, and D = 2 scheme choice. 



Figure 4. \V US \ vs. so for a selection of AHm 
FESRs and, for comparison, the CIPT+Adler w T 
AIE determination. 



4. APPENDIX: Numerical results for the 



OPE input to the FB s 



FESR 



This appendix contains numerical results for 
the OPE contributions to the conventional inclu- 
sive, non-longitudinally-subtractcd, FB, sq = 
FESR determination of \V US \, in a form that al- 
lows the reader to perform further explorations. 
For this case, in the notation of Ref. [5] , one has 



\V US 



Rt,V + A 

\v ud \ 2 



SRr. 



(10) 



th 



Table 1 

8Rope d=2 ^ 0r var ious tit and D = 2 evalua- 
tion schemes (CI+Ad, CI+co and FOPT denote 
the CIPT+Adler function, CIPT+correlator and 
FOPT schemes, respectively) 



Scheme 



Tlx 



SR 



L+T 

OPE:D=2 



CI+Ad 


3 


0.0578 + 0.0013 




4 


0.0494 + 0.0025 


CI+co 


3 


0.0515 + 0.0022 




4 


0.0371 + 0.0040 


FOPT 


3 


0.1002 + 0.0016 




4 


0.1046 + 0.0019 



With ChPT input for m s /(m d +m u ), GMOR 
for {{md+m u )uu} , and the input already specified 
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above, the D = 4 contribution is given by 
6r op T e-d=4 = 0.0072(2)r c - 0.00030 

with r c = (ss)/(uu), and the three terms on the 
RHS being associated with the strange quark con- 
densate, light quark condensate and 0(mf ) con- 
tributions, respectively. The errors shown are 
those associated with the uncertainty in the in- 
put a s . The value r c = 0.8 employed in earlier 
analyses, follows from a 2002 sum rule analysis 
of /b 3 I Ib with then-current quenched lattice re- 
sults for the decay constant ratio as input [2~7] . 
Current n.j = 2+1 results for this ratio yield in- 
stead r c = 1.2. Finally, if desired, an estimate 
of the small D = 6 contribution can be included 
using a rescaled version of the vaccum saturation 
approximation (VSA). This yields 

6RoP%D- e = -0.0043(l-r c 2 ) ( . ) 
OPE D - e v c > \ 0.000150 GeV 6 / 

(13) 

where p parametrizes VSA breaking. The value in 
the denominator of the last factor reflects Ioffe's 
assessment of p [26] . 
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